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Abstract. A novel approach to electronic correlations in magnetic crystals which takes 
into account a dynamical many-body effects is present. In order to to find a frequency 
dependence of the electron self energy, an effective quantum-impurity many-particle 
problem need to be solved within the dynamical mean-field theory. The numerically 
exact QMC-scheme and the spin-polarized fluctuation exchange approximation are 
used for the self-consistent solution of this single-site many-particle problem. The cal- 
culations of effective exchange interaction parameters based on the realistic electronic 
structure of correlated magnetic crystals have been discussed. 



1 Introduction 

The calculation of thermodynamic properties and excitation spectra of different 
magnetic materials is one of the most important problems of the microscopic 
theory of magnetism. A modern computational material science is based on the 
density functional (DF) approach |Q. It is a common practice to use this scheme 
not only for the total energy calculations and related quantities such as charge 
and spin densities, but also for different spectral characteristics. Sometimes an 
agreement of the computational results with the experimental data is very im- 
pressive, despite the absence of a reliable theoretical background. In principle, 
the energies of Kohn-Sham quasiparticlies Q which are calculated in a stan- 
dard band theory are just auxiliary quantities for the total energy calculation. 
Among all achievements of the quantitative electronic theory a list of difficulties 
and shortcomings grows, especially when one consider the magnetic d- and f- 
electron systems. In a number of cases the theory appeared to be qualitatively 
inadequate. First, the DF scheme cannot describe correctly the phenomenon of 
"Mott insulators" as it was first shown by Terakura et al Q in attempt to 
calculate the electronic structure of 3d-metal oxides. Later we faced similar prob- 
lems in field of the high-Tc superconductors [Q and other compounds § . The 
Ce- and U- based "heavy fermion" compounds such as CeCug, UPta, etc, are an- 
other "hot-spots" : normally the calculated values of effective masses are orders 
of magnitude smaller than the experimentally observed ones Even for the 
pure 3d metals some qualitative differences between the theory and experiments 
exists. For example, there are at least three difficulties with the photoelectron 
spectra of ferromagnetic nickel : (i) the measured width of the occupied part 
of d-band is 30% narrower than the calculated one (ii) the spin splittings is 
twice smaller than in the LSDA and (iii) the band structure cannot describe 
a famous 6 eV - satellite. Calculations for paramagnetic spin-disordered states 
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[|j lead to the conclusion that Ni has no local magnetic moments above the 
Curie temperature Tc, in a clear contradiction with experimental results For 
iron, a standard band theory cannot explain the data about spin polarization of 
thermionic electrons [ p^PTl| , p^ and some features of angle-resolved photoemision 
spectra [ p3| , p"4p^ . All these difficulties put many questions to the DF approach: 
what are the "electron spectrum" which we really calculated and how one can 
improve the electronic theory for magnetic d- and f- systems? 

It was understood many years ago that all this problems are connected with 
inadequate description of many-body effects in DF calculations of the excita- 
tion spectra. Such methods as GW (l^, LDA-I-U ^ and SIC (Self- Interaction 
Correction approximation) |16 have been proposed to improve the situation. 
This methods are very useful for the description of antiferromagnetic transition- 
metal oxides as the Mott insulators (l^,^. However, one should note that both 
LDA-I-U and SIC are just the mean-field approximations and cannot describe 
the correlation effects which are, by definition, the many-body effects beyond 
the Hartree-Fock. For example, in these approaches one need spin or orbital or- 
dering to describe the Mott insulator and it is impossible to describe correctly 
electronic structure of NiO or MnO in paramagnetic phase. At the same time, 
the magnetic ordering should not be important for the basic physics of Mott 
insulators |^ . All the "Hartree-Fock- like" approaches fail to describe the renor- 
malization of the effective mass in the heavy fermion systems. There are also 
many problems concerning the electronic structure and itinerant magnetism of 
3d metals as described above. Thus, one need some practical ways to incorporate 
correlation effects in the electronic structure of solids. 

In principle, there are two ways to include them into DF calculations. The 
first one is the use of time- dependent DF formalism which can guarantee, in prin- 
ciple, an opportunity to calculate exact response functions l]!^ , in the same sense 
as the Hohenberg-Kohn theorem guarantees the total energy in usual "static" 
DF However, all the expressions for this time-dependent non-local DF in 
real calculations are based on RPA-like approximations which described not 
satisfactory the really highly correlated systems. They are excellent for investi- 
gation the plasmon spectrum of aluminum, but not for understanding a nature 
of high- Tc superconductivity or the heavy fermion behavior. Another way is to 
use an "alternative" many-body theory developed in the 50-th by Gell-Mann 
and Brueckner, Galitskii and Migdal, Beliaev and many others in terms of the 
Green functions rather than in the electron density jl^. We try to formulate 
such computational approach as a generalization of LDA-f U scheme, a so-called 
"LDA-I--I-" method. The main difference of LDA++ approach from LDA+U one 
is an account in the former dynamical fiuctuations, or the real correlation ef- 
fects, described by local but energy dependent self-energy S{E)^ so the LDA-|--|- 
means LDA + U + S{E). 

The comparison of the standard DF theory in the local spin density approx- 
imation (LSDA) and LDA-I--I- approach is represented in the table I. 

First of all the LSDA theory is based on the Hohenberg - Kohn theorem 
that the total energy Etot is a functional of charge and spin densities, while the 
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Table 1. Comparison of different schemes 



LSDA 



LDA++ 



Density functional 
Density p(r) 
Potential Vxcir) 



Etot 



Edc 



Edc = Eh + J pVxcdr ~ E^c 



Baym-Kadanoff functional 
Green- Function G(r, r', E) 
Self-energy E,{E) 

lT? — ^ sp ~ ^dc 

Osp = -Trln[-G-^] 
Ode = TrSG ~ <Plw 



LDA+-f scheme considers the thermodynamic potential Q as a. functional of ex- 
act one-particle Green functions. This approach in many-particle theory has been 
introduced in the works by Luttinger and Ward and Baym and Kadanoff 
[ pO| . The Green function in LDA-I--I- theory plays the same role as the density 
matrix in LSDF formalism. We stress the dynamic nature of the correlation ef- 
fects which are taken into account in the LDA-|— I- approach since the density in 
the LSDA is just the static limit of the local Green function. Further, the self 
energy S is analogous to the exchange-correlation potential; local approximation 
for S, which is assumed to be energy-dependent but not momentum-dependent 
corresponds to the local approximation for V^c- In both formalisms the thermo- 
dynamic potential can be represented as a "single-particle" one, H^p minus the 
contributions of the "double counted terms", f^dc- It will be important for the 
consideration of so-called "local force theorem" and the computation of mag- 
netic interaction parameters (Sect. III). The single-particle contribution to the 
thermodynamic potential in the LDA-I--I- would have the same form as in the 
LSDA if we take into account only a pole part of the Green function and neglect 
the quasiparticle damping. However, even in that case the quasiparticle ener- 
gies are not the same since the poles of the Green functions are not coincide, 
generally speaking, with the "Kohn-Sham" energies. The quantity ^lw is the 
Luttinger- Ward generating functional for the self energy, or the sum of all the 
skeleton diagrams without free legs ]l9t . 

The difficulty with the finite temperature effects is one of the main shortcom- 
ings of a standard DP formalism. In all realistic calculations the temperature is 
included in the Fermi distribution functions and in the lattice constants via the 
thermal expansion [ pTf . At the same time, for the itinerant electron magnets the 
temperature effects connected with the "Bose" degrees of freedom due to spin 
waves and paramagnons are much more important ||2^ . In principle, these effects 
could be taken into account in the DP theory via the temperature dependence 
of the exchange-correlation potential, the corresponding terms being nonlocal. It 
is not easy to propose an adequate expression for such temperature-dependent 
non-local potential. One of the first attempts in this direction based on simple 
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RPA-like considerations [g3|. On the other hand, in LDA++ type of scheme all 
the calculations are naturally carried out for finite temperatures by the using of 
Matsubara frequencies, as in the usual many-body theory Q. 

The main assumption of the LDA-I--I- approach is the importance of only in- 
trasite "Hubbard correlations" with the local approximation for the self-energy. 
It is worthwhile to stress a difference of this kind of locality from the locality 
in DF theory. In the latter, the local approximation means that the exchange- 
correlation energy is calculated for the homogeneous electron gas . It is known 
from exact QMC calculations that the correlation effects could lead to some in- 
stabilities of the state of homogeneous electron gas (magnetism, charge ordering, 
etc) only for electron densities which are order of magnitude smaller than ones 
typical for real metals (the critical values of the parameter rg are of order of hun- 
dred in comparison with the "normal" range 2-6 for metals). At the same time, 
magnetism and charge ordering are rather usual for real compounds with the 
d- and f- elements. It seems that the "atomic-like" features of d- and f- states 
are of the crucial importance to describe the correlation effects in real com- 
pounds. Only these features are taken into account in the Hubbard-like terms 
for the d- or f-states in LDA-I--I- approach. Therefore one can view the LDA-I--I- 
as a simplest way for quantitative considerations of the correlation effects in the 
transition metals, and their compounds, based on the LSDA description for all 
non-correlated electrons in the systems. 

The investigation of correlation effects in the electronic structure and mag- 
netism of iron-group metals is still far from the final picture and attracts contin- 
uous interest (see, e.g., [p|, p5|j2^ , |27[ | and Refs therein). Despite of many attempts, 
the situation is still unclear both theoretically and experimentally. For example, 
there is no agreement on the presence of 5 eV satellite in photoemission spectrum 
of iron , and on the existence of local spin splitting above Curie tempera- 

ture in nickel . From the theoretical point of view, different approaches such 
as the second-order perturbation theory [ p9|j26[, the T-matrix approximation 
[ psIJsoll , the three-body Faddeev approximation^^, and the moment expansion 
method |^ were used. Unfortunately, the applicability of these schemes are 
not clear. Here we present the version of "LDA-f-f" approach [|2|j3^,Q which 
is based on the combination of standard band theory technique with so-called 
dynamical mean-field theory (DMFT) or LDA+DMFT scheme fUl . 

2 Computational technique 

We start from LDA+U Hamiltonian in the diagonal density approximation: 
- V t^^^ ,c+ cv , +i V IP ,n n- , + 



2 5Z ^mm' J mm'^'^ima'^im' a (1) 
im^m' a 

where i is the site index and m is the orbital quantum numbers; cr — I, | is the spin 
projection; c+,c are the Fermi creation and annihilation operators [n — c^c); 
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t is effective single-particle Hamiltonian obtained from the non-magnetic 
LDA with the correction for double counting of the average interactions among d- 
electrons. In the general case of spin-polarized LSDA Hamiltonian this correction 
is presented in Refs. ||3^,|l^,^. In the non-magnetic LDA this is just a shift 
"back" of correlated d-states with respect to s,p-states by the average Coulomb 
and exchange potential: Ad = U(nd ~ ^) ~ ■^Jind — 1), where U and J are the 
average values of Umm' and Jmm' matrices and Ud is the average number of 
d-electrons. 

The screened Coulomb and exchange vertex for the d-electrons 

Umm' =< mm'|y//^(r - r')\mm' >, Jw =< mm'\Vf^l{r - r')\m'm > (2) 

are expressed via the effective Slater integrals. We use the minimal spd-basis 
in the LMTO-TB formalism and numerical orthogonalization for i-^^^(k) 
matrix |33j. Local density approximation was used for the self-consistent 
electronic structure calculation. 

In order to find the best local approximation for the self-energy we use the 
DMFT method |Q for real systems. This scheme become exact in the limit of 
infinite lattice coordination number |^ . The DMFT approach reduce the lattice 
many-body problem (Eq.(||)) to the self-consistent solution of effective one-cite 
Andersen model. In this case we need a local Green-function matrix which has 
the following form in the orthogonal LMTO-representation: 

G{iLu) = ^{zt^ + fi- t^°^{k) - S{iu;)}-^ (3) 

k 

were fi is the chemical potential. Note that due to cubic crystal symmetry of 
ferromagnetic bcc-iron the local Green function without spin-orbital interactions 
is diagonal both in the orbital and the spin indices. The so-called bath Green 
function which defined the effective Andersen model and preserve the double- 
counting of the local self-energy is obtained as a solution of the following impurity 
model ||: 

Go„\M = G,-iM + ^™M (4) 

The local Green functions for the imaginary time interval [0, /3] with the 
mesh Ti = I At, I = 0, L — 1, and At ~ (3/L, where /3 = ^ is calculated in the 
path- integral formalism |3|,^ : 

G!i' = |^det[0(s)]*G^:(.) (5) 

here we redefined for simplicity m = {m,a},Z is the partition function and 
the so-called fermion-determinant det [0(s)] as well as the Green function for 
arbitrary set of the auxihary fields G(s) = O^^(s) are obtained via the Dyson 
equation |^ for imaginary-time matrix (G„i(s) = G''J^{s)): 

G,„ = [l-(GO,-l)(e^"-iriG°, 
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where the effective fluctuation potential from the Ising fields sl^^, = ±1 is 



T/' — \ ' ti _j- l,m<m' 

m. ~^ / , '^mm' Sjyi^i CTmm' , WnCre (Jmm' — i m > m! 

and the discrete Hubbard-Stratonovich parameters are A„i„i' — arccosh[exp(^Z\Tt/„ 

Using the output local Green function from QMC and input bath Green 
functions the new self-energy is obtain via Eq.(Q) and the self-consistent loop 
can be closed through Eq.(||). The main problem of the multiband QMC for- 
malism is the large number of the auxiliary fields sj„„j, . For each time slices I 
it is equals to M{2M — 1) where M is the total number of the orbitals which 
gives 45 Ising fields for the d-states case. We compute the sum over this auxil- 
iary fields in Eq. (^) using important sampling QMC algorithm and performed a 
dozen of self-consistent iterations over the self-energy Eqs.(|[||j|). The number of 
QMC sweeps was of the order of 10^ on the CRAY-T3e supercomputer. The final 
G„i(t) has very little statistical noise. We use maximum entropy method for 
analytical continuations of the QMC Green functions to the real axis. Compar- 
ison of the total density of states (DOS) with the results of LSDA calculations 
(Fig.^ shows a reasonable agreement for single-particle properties of not "highly 
correlated" ferromagnetic iron. We calculate the bcc iron at experimental lattice 
constant with 256 k-points in the irreducible part of Brilloiun zone. The Mat- 
subara frequencies summation corresponds to the temperature of about T=850 
K. The average magnetic moment is about 1.9 /is which corresponds to a small 
reduction of the LSDA- value of 2.2 fiB for such a high temperature. The DOS 
curves in the LDA-I-Z" approach with exact QMC solution of on-site multiorbital 
problem is similar to that obtained within the simple perturbative fluctuation- 
exchange (FLEX) approximation described below. The discussion of the results 
and the comparison with the experimental data will be given in Section 4. 

The QMC method described above is probably the most rigorous real way to 
solve an effective impurity problem in the framework of DMFT theory. However, 
it is rather time consuming. Besides that, in the previous section we did not work 
with complete four-indices Coulomb matrix: 



(12|t;|34) = J drdr'riir)r2ir')v,erir-r')Mr)Mr'), (6) 

where we define for simplicity mi = 1. 

For moderately strong correlations (which is the case of iron group metal) 
one can propose an approximate scheme which is more suitable for the calcula- 
tions. It is based on the fluctuation exchange (FLEX) approximation by Bickers 
and Scalapino generalized to multiband spin-polarized case p3| , p^p^ . The 
electronic self-energy in the FLEX is equal to: 



(7) 
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Fig. 1. Spin-resolved density of d-states and magnetic moments for ferromagnetic iron 
in the LSDA and the LDA+QMC calculations for different average Coulomb interac- 
tions with J=0.9 eV and temperature T=1500 K. 

where the Hartree-Fock contribution has a standard form: 



y^HF 

^12, (T 



E 

34 



(13|H24)^n-;-(13|H42X4 



(8) 



with the occupation matrix ~ G'2i(t — > —0); this contribution to E is 
equivalent to spin-polarized "rotationally-invariant" LDA-I-U method p6[ . 
The second-order contribution in the spin-polarized case reads: 



.(2) 



W = - E (13 74) G?8(r)* 

{3-8} 



(85 \v\ 26) J2 G^3MGI;(-t) - (85 \v\ 62) G^,,{T)Gl,i-T) 
and the higher-order particle-hole (or particle-particle) contribution 

34, cr' 



with p-h (p-p) fluctuation potential matrix: 



H/TT y/U (it^) 
H/J-f (iw) W^^ {iiJ) 



(9) 



(10) 



(11) 
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where the spin-dependent effective potentials has a generahzed RPA-form and 
can be found in Note that for both p-h and p-p channels the effective 
interactions, according to Eq.(|l^), are non-diagonal matrices in spin space as 
in the QMC-scheme, in contrast with any mean-field approximation like LSDA. 
This can be important for the spin-dependent transport phenomena in transition 
metal multilayers. 

We could further reduce the computational procedure by neglecting dynam- 
ical interaction in the p-p channel since the most important fluctuations in itin- 
erant electron magnets are spin-fluctuations in the p-h channel. We take into 
account static ( of T- matrix type) renormalization of effective interactions re- 
placing the bare matrix 1/12,34 = (12 \v\ 34) in FLEX-equations with the corre- 
sponding spin-dependent scattering T— matrix 



12 



rj-iO 



34) = (12 \v\ 34) - ^ (12 \v\ 56) / drG^g (r) G?; (r) (78 

5678 



34 ~ 



(12) 



Similar approximation has been checked for the Hubbard model pq | and 
appeared to be accurate enough for not too large U. Finally, in the spirit of 
DMFT-approach U = U[Go], and all the Green functions in the self-consistent 
FLEX-equations arc in fact the bath Green- functions Gq. 



3 Exchange interactions 

An useful scheme for analyses of exchange interactions in the LSDF approach 
is a so called "local force theorem". In this case the calculation of small total 
energy change reduces to variations of the one-particle density of states p^ , p7[ |. 
First of all, let us prove the analog of local force theorem in the LDA+-I- ap- 
proach. In contrast with the standard density functional theory, it deals with the 
real dynamical quasiparticles defined via Green functions for the correlated elec- 
trons rather than with Kohn-Sham "quasiparticles" which are, strictly speaking, 
only auxiliary states for the total energy calculations. Therefore, instead of the 
working with the thermodynamic potential as a density functional we have to 
start from the general expression for f2 in terms of exact Green function in the 
Table I. We have to keep in mind also Dyson equation = Gq^ — S and 
the variational identity S^Plw = TrSSG. Here Tr = Tr^iL„ is the sum over 
Matsubara frequencies Tr^... — TJ2 ^ — t^T {2n -t- 1) , n = 0, ±1, T is the 

temperature, and iLa are site numbers (i), orbital quantum numbers (L = Z,rn) 
and spin projections a , correspondingly. We represent the expression for ]7 as 
a difference of "single particle" (sp) and "double counted" (dc) terms as it is 
usual in the density functional theory. When neglecting the quasiparticle damp- 
ing, J7sp will be nothing but the thermodynamic potential of "free" fermions 
but with exact quasiparticle energies. Suppose we change the external potential, 
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for example, by small spin rotations. Then the variation of the thermodynamic 
potential can be written as 

5Q = S* + SiQsp — 5 fide (13) 

where 5* is the variation without taking into account the change of the "self- 
consistent potential" (i.e. self energy) and 5i is the variation due to this change 
of S . To avoid a possible misunderstanding, note that we consider the variation 
of fi in the general "non-equilibrium" case when the torques acting on spins are 
nonzero and therefore 5f2 ^ 0. In order to study the response of the system to 
general spin rotations one can consider either variations of the spin directions 
at the fixed effective fields or, vice versa, rotations of the effective fields, i.e. 
variations of at the fixed magnetic moments. We use the second way. Taking 
into account the variational property of one can be easily shown (cf. Ref. p^Oj) 
that SiQsp = 5^dc = TrGSS and hence 

5Q = 5*^2sp = -5*Tr\n[lJ -Gq^] (14) 

which is an analog of the " local force theorem" in the density functional theory 

Further considerations are similar to the corresponding ones in LSDF ap- 
proach. In the LDA-I--I- scheme, the self energy is local, i.e. is diagonal in site 
indices. Let us write the spin-matrix structure of the self energy and Green 
function in the following form 

Si = Si + Sla , Gij = G^j + G -jO- (15) 

where E^^'^'"^ ~ ^ (^sj ± S^^, Sf — ^fe^, with e.^ being the unit vector in the 

direction of effective spin-dependent potential on site i,a = {a^, ctj,, cr^) are Pauli 
matrices, G^^ = ^Tr„{Gij) and G*^ = jTr„{Gija). We assume that the bare 
Green function does not depend on spin directions and all the spin-dependent 
terms including the Hartree-Fock terms are incorporated in the self energy. Spin 
excitations with low energies are connected with the rotations of vectors e^: 
5ei = 5iPi X Bi According to the "local force theorem" (|lj) the corresponding 
variation of the thermodynamic potential can be written as 5Q — 5* Qgp = 'ViSipi 
where the torque is equal to 

V,; = 2Tr^L [St X G^J (16) 

Using the spinor structure of the Dyson equation one can write the Green 
function in this expression in terms of pair contributions (a similar trick has 
been proposed in Ref. |Q in the framework of LSDF approach). As a result, 
we represent the total thermodynamic potential of spin rotations or the effective 
Hamiltonian in the form ||3^ 

i^spin = — ""^Tr^L {{^ijSj) i^ji^i) ~ SiGijSjG'ji — i (S^ x G^jS^j) Gjj} 

(17) 
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Vi This means that 

G— const 



one can show by direct calculations that 

^spin {^i} is the effective spin Hamiltonian. The last term in Eq.([l7[) is nothing 
but Dzialoshinskii- Moriya interaction term. It is non-zero only in relativistic case 
where Sj and GJj can be, generally speaking, "non-parallel" and dj ^ Gji for 
the crystals without inversion center. 

In the nonrelativistic case one can rewrite the spin Hamiltonian for small 
spin deviations near coUincar magnetic structures in the following form fispin = 
~ Jij^i^j where 

J^J = -Tr^L {s!Glsp\^ (18) 

are the effective exchange parameters. This formula generalize the LSDA expres- 
sions of to the case of correlated systems. 

Spin wave spectrum in ferromagnets can be considered both directly from 
the exchange parameters or by the consideration of the energy of corresponding 
spiral structure (cf. Ref. ^^). In nonrelativistic case when the anisotropy is 
absent one has 

^ ^ E ^OJ- (1 - cosqR^.) ^ ^[J(0) - J(q)] (19) 

j 

where M is the magnetic moment (in Bohr magnetons) per magnetic ion. 

It should be noted that the expression for spin stiffness tensor Daf3 defined by 
the relation ujq = Dapqaqp (q 0) in terms of exchange parameters has to be 
exact as the consequence of phenomenological Landau- Lifshitz equations which 
are definitely correct in the long- wavelength limit. Direct calculation basing on 
variation of the total energy under spiral spin rotations (cf. Ref. |Q) leads to 
the following expression 

were k is the quasimomentum and the summation is over the Brillouin zone. The 
expressions Eqs.(p^) and ( |T9| ) are reminiscent of usual RKKY indirect exchange 
interactions in the s-d exchange model (with U'^ instead of the s-d exchange 
integral) . 

We prove in the Appendix that the expression for the stiffness is exact within 
the local approximation. At the same time, the exchange parameters themselves, 
generally speaking, differ from the exact response characteristics defined via 
static susceptibility since the latter contains vertex corrections. The derivation 
of approximate exchange parameters from the variations of thermodynamic po- 
tential can be useful for the estimation of in the different magnetic systems. 



4 Computational results 



We have started from the spin-polarized LSDA band structure of ferromagnetic 
iron within the TB-LMTO method in the minimal s,p, d basis set and used 
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numerical orthogonalization to find the Ht part of our starting Hamiltonian. 
We take into account Coulomb interactions only between d-states. The correct 
parameterization of the Hjj part is indeed a serious problem. For example, first- 
principle estimations of average Coulomb interactions (U) in iron lead 
to unreasonably large value of order of 5-6 eV in comparison with experimental 
values of the U-parameter in the range of 1-2 eV Semiempirical analysis 
of the appropriate interaction value Q gives U ~ 2.3 eV. The difficulties with 
choosing the correct value of U are connected with complicated screening prob- 
lems, definitions of orthogonal orbitals in the crystal, and contributions of the 
intersite interactions. In the quasiatomic (spherical) approximation the full U- 
matrix for the d— shell is determined by the three parameters U, J and 5 J or 
equivalently by effective Slater integrals F°, and F^ ||3^,P- For example, U= 
F°, J=(F^+F^)/14 and we use the simplest way of estimating &J or F^ keeping 
the ratio F^/F^ equal to its atomic value 0.625 |5l| . 

Note that the value of intra-atomic (Hund) exchange interaction J is not 
sensitive to the screening and approximately equals to 0.9 eV in different es- 
timations For the most important parameter U, which defines the bare 
vertex matrix Eq.(||), we use the value U = 2.3 eV for Fe j50|, C/ = 3 eV for 
Co and Mn and [7 = 4 eV for Ni and Cu. To calculate the spectral functions: 
Acr (k, E) — —^Tr^Ga (k, E + iO) and DOS as their sum over the Brillouin zone 
we first made analytical continuation for the matrix self-energy from Matsub- 
ara frequencies to the real axis using the Fade approximation [^2[ , and then 
numerically inverted the Green-function matrix as in Eq. (||) for each k-point. 
In the self-consistent solution of the FLEX equations we used 1024 Matsubara 
frequencies and the FFT-scheme with the energy cut-off at 100 eV. The sum over 
irreducible Brillouin zone have been made with 72 k-points for SCF-iterations 
and with 1661 k-points for the final total density of states. 

The depolarization of states near the Fermi level is another important correla- 
tion effect. The decrease of the ratio P = [NiIEf) - (Ep)] / [N^ {Ep) + A^x {Ep)] 
is a typical sign of spin-polaron effects [P7y53[| . In our approach this effects 
are taken into account through the W^^^^^ terms in the effective spin-polarized 
LDA-h-h potential (Eq. (0)). 

The energy dependence of self-energy in Fig.^ shows characteristic features of 
moderately correlated systems. At low energies \E\ < 1 eV we see a typical Fermi- 
Hquid behavior ImU (E) —E^, dReS (E) /dE < 0. At the same time, for the 
states beyond this interval within the d-bands the damping is rather large (of the 
order of 1 eV) so these states corresponds to ill-defined quasiparticles, especially 
for occupied states. This is probably one of the most important conclusions of 
our calculations. Qualitatively it was already pointed out in Ref. |^ on the basis 
of a model second-order perturbation theory calculations. We have shown that 
this is the case of realistic quasiparticle structure of iron with the reasonable 
value of Coulomb interaction parameter. 

Due to noticeable broadening of quasiparticle states the description of the 
computational results in terms of effective band structure (determined, for ex- 
ample, from the maximum of spectral density) would be incomplete. We present 
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Fig. 2. Total spin-polarized density of states and d-part of self-energy for iron with 
U=2.3 eV and J=0.9 eV for the temperature T=750 K. Two different self-energies for 
Xig and eg d-states in the cubic crystal field symmetry are presented and four different 
lines corresponds to imaginary part spin-up (full line) and spin-down (dashed line) as 
well as real part spin-up (dashed-dot line) and spin-down (dashed-double-dot line). 



on the FigJ^ the full spectral density Aa (k, E) including both coherent and inco- 
herent parts as a function of k and E. We see that in general the maxima of the 
spectral density (dark regions) coincide with the experimentaUy obtained band 
structure. However, for occupied majority spin states at about -3 eV the distri- 
bution of the spectral density is rather broad and the description of this states 
in terms of the quasiparticle dispersion is problematic. This conclusion is in 
complete quantitative agreement with raw experimental data on angle-resolved 
spin-polarized photoemission |Q with the broad non-dispersive second peak in 
the spin-up spectral function around -3 eV. 

Comparison of the DOS for transition metals in the Fig. ^ shows an inter- 
esting correlations effects. First of all, the most prominent difference from the 
LSDA calculation is observed for the antiferromagnetic fcc-Mn. There is clear 
formation of the lower and upper Hubbard bands around ± 3 eV. Such behavior 
is related with the half-field Mn d-shell which corresponds to a large phase space 
for particle-hole fluctuations. For the ferromagnetic bcc-Fe the p-h excitations 
are suppressed by the large exchange splitting and a bcc structural minimum in 
the DOS near the Fermi level. In the case of ferromagnetic fcc-Co and Ni the 
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Fig. 3. Spectral function of ferromagnetic iron for spin-up (a) and spin-down (b) and 
the two k-directions in the Brillouin zone compare with the experimental angle resolved 
photoemission and de Haas - van Alphen (at the Ei?=0) points 



correlation effects are more important then for Fe, since there is no structural 
bcc-dip in the density of states. One could see a formation of a "three-peak" 
structure for the spin-down DOS for Co and Ni and satellite formation around 
-5 eV. In order to describe the satellite formation more carefully one need to 
include T-matrix effects or use the QMC scheme in LDA-I-DMFT calcu- 

lations. Finally, there is no big correlation effects in non-magnetic fcc-Cu, since 
the d-states are located well bellow the Fermi level. 

Using the self-consistent values for computed by QMC technique 

(Section II) we calculate the exchange interactions (Eq.|l^ and spin-wave spec- 
trum (Eq.^|) using the four-dimensional fast Fourier transform (FFT) method 
for (k, iuj) space with the mesh 20"^ x 320. The spin-wave spectrum for 
ferromagnetic iron is presented in Fig.^ in comparison with the results of LSDA- 
exchange calculations 1 47 and with different experimental data |^ , ^7|j5^ . This 
room-temperature neutron scattering experiments has a sample dependence (Fe- 
12%Si in Ref. and Fe-4%Si in Ref. Hi) due to problems with the bcc-Fe 

crystal growth. Note that for high-energy spin-waves the experimental data |Q 
has large error-bars due to Stoner damping (we show one experimental point 
with the uncertainties in the q space). On the other hand, the expression of 
magnon frequency in terms of exchange parameters itself becomes problematic 
in that region due to breakdown of adiabatic approximation, as it is discussed 
above. Therefore we think that comparison of theoretical results with experi- 
mental spin- wave spectrum for the large energy needs additional investigation of 
Stoner excitation and required calculations of dynamical susceptibility in the 
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Fig. 4. The spin-polarized partial 3d density of states for different transition metals 
at the temperature T=750 K. The full line is the spin-up DOS and the dashed line is 
the spin-down DOS. 



LDA-f-f approach p8[. Within the LSDA scheme one could use the linear- 
response formalism [ ^9| to calculate the spin-wave spectrum with the Stoner 
renormalizations, which should gives in principle the same spin- wave stiffness 
as our LSDA calculations. Our LSDA spin-wave spectrum agree well with the 
results of frozen magnon calculations l6Q,pl| . 



At the lower-energy, where the present adiabatic theory is reliable, the LDA++ 
spin-waves spectrum agree better with the experiments then the result of the 
LSDA calculations. Experimental value of the spin-wave stiffness D=:280 meV/A^ 
agrees well with the theoretical LDA++ estimations of 260 meV/A^ ]3^ . 
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Fig. 5. The spin- wave spectrum for ferromagnetic iron in the LSD A and LDA+X" 
approximations compared with different experiments (circles [16], squares [17], and 
diamonds [18]) (a); The corresponding spin- wave spectrum from LDA-I-I7 scheme in 
the (110) plane (b). 

Appendix 

Here we prove that the expression for the stiffness constant Eq.(^) is exact in 
the framework of DMFT scheme. 

A rigorous expression for the stiffness constant has been obtained in Ref. 



2<S'' > D^pq^qp = C{(i)-i<l- / ^(fc A: t)Gx(fc)GT(fc)(vk • q) (21) 



k = (k, fco) is a momentum-energy 4- vector, G{k) is the electron Green function, 
t(k) is the one-particle (band) Hamiltonian, n(k) is the one-particle distribution 
function, 7 is the irreducible vector three-leg vertex. The irreducible scalar and 
vector vertices, 7 and 7, are connected with the total ones, Fq and F, by Bethe- 
Salpeter equations and the total vertices satisfy the Ward-Takahashi identity 
which is the consequence of the rotational invariance. Taking into account the 




where 




16 A. I. Lichtenstein et al. 



Bethe-Salpeter equations one can rewrite them in terms of the irreducible ver- 



tices (for more details, see |62 ). Let us assume now that the scalar irreducible 
vertex 7 and S are momentum-independent (which is the real case in DMFT 
approach). Considering the q-dependent part of the Ward-Takahashi identity in 
the limit q — > we have: 7(fc i,k ]) — = v(k). Substituting this into 



Eq.(|2lD one has 



On the other hand, our Eq.(f9|) reads 2 < S"^ > Dap = k dq'df^ U-^Q where 



J(q) is defined by Eq.(|l8|). Calculating the derivatives of the exchange parame- 
ters we obtain 

(we shift here k — > k + q in the integrand) . Then 

This expression can be simplified by using the sum rule (G^ — G j ) = G| SG^^ 
(where S = ^ ^l) which is the consequence of the Dyson equation provided 



that only S is spin-dependent. Taking it into account in Eq.(22) one has 



'''''' -7(0-''-'<^'1?-"^'l?-i:-"<''t^t'' 



c^ag^g " J (27r)4 ' dka, ^ dkp ' dka dkfi 

The first term is exactly coincide with the last one in Eq.(^l]) and the first 
•m can be transformed further using the identity: Gl{k) = Then 

— 9t(k)9t(k) ^ 9n^(k) at(k) ai(k) 



[ d^k ^ , 9t(k)9t(k) _^ 



(27r)4i ^ ' dka dkp ^ at(k) dk^ dkp 

Since nCT(k) =ncr\t{}i)\ we have: g|-no.[i(k)] = ^dt(^) %k^'' and finally, integrating 
by part, one obtains: 



dn„ dt(k) ^ ^ dH{k) 



_ V^^l^ - Vn (k) 

^ dka dkp ^ 'dk^dkf} 

Thus, our expression (Eq.([20|)) coincides with the exact one (Eq.(^l|)). We use 
here the only assumption that both the self-energy and three-leg irreducible 
vertex are momentum independent as well as the Ward-Takahshi identities which 
are exact consequences of the rotationally invariance of the spin system. 
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5 Conclusions 

We have proposed a general scheme for investigation of the correlation effects in 
the quasiparticle band structure calculations for itinerant-electron magnets. This 
approach is based on the combination of the dynamical mean-field theory and 
the fluctuating exchange approximation. Application of LDA+DMFT method 
gives an adequate description of the quasiparticle electronic structure for fer- 
romagnetic iron. The main correlation effects in the electron energy spectrum 
are strong damping of the occupied states below 1 eV from the Fermi level Ep 
and essential depolarization of the states in the vicinity of Ep. We obtained a 
reasonable agreement with different experimental spectral data (spin-polarized 
photo- and thermoemission). The method is rather universal and can be applied 
for other magnetic systems, both ferro- and antiferromagnets. 

We discussed as well a general method for the investigation of magnetic in- 
teractions in the correlated electron systems. This scheme is not based on the 
perturbation theory in "[/" and could be applied for rare-earth systems where 
both the effect of the band structure and the multiplet effects are crucial for a 
rich magnetic phase diagram. Our general expressions are valid in relativistic 
case and can be used for the calculation of both exchange and Dzialoshinskii- 
Moriya interactions, and magnetic anisotropy jsj]. An illustrative example of 
ferromagnetic iron shows that the correlation effects in exchange interactions 
may be noticeable even in such moderately correlated systems. For rare-earth 
metals and their compounds, colossal magnetoresistance materials or high- Tc 
systems, this effect may be much more important. For example, the careful in- 
vestigations of exchange interactions in MnO within the LSDA, LDA-t-U and 
optimized potential methods for MnO show the disagreement with exper- 
imental spin- wave spectrum (even for small q), and indicate a possible role of 
correlation effects. 

This work demonstrates an essential difference between spin density func- 
tional approach and LDA++ formalism. The latter method deals with the ther- 
modynamic potential as a functional of the local Green function rather than the 
electron density. Nevertheless, there is a close connection between two techniques 
(the self-energy corresponds to the exchange- correlation potential, etc) . In par- 
ticular, an analog of local force theorem can be proved for LDA-|--|- approach. It 
may be useful not only for the calculation of magnetic interactions but also for 
elastic stresses, in particular, pressure, and other physical properties. 
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